Draft version November 10, 2008 

Preprint typeset using I£T^X style cmulateapj v. 10/09/06 



IONIZING PHOTON ESCAPE FRACTIONS FROM HIGH REDSHIFT DWARF GALAXIES 

John H. Wise 1,2 and Renyue Cen 3 

Draft version November 10, 2008 

ABSTRACT 

It has been argued that low-luminosity dwarf galaxies are the dominant source of ionizing radiation 
during cosmological reionization. The fraction of ionizing radiation that escapes into the intergalactic 
medium from dwarf galaxies with masses less than ~10 95 solar masses plays a critical role during 
this epoch. Using an extensive suite of very high resolution (0.1 pc), adaptive mesh refinement, 
radiation hydrodynamical simulations of idealized and cosmological dwarf galaxies, we characterize 
the behavior of the escape fraction in galaxies between 3 x 10 6 and 3 x 10 9 solar masses with different 
spin parameters, amounts of turbulence, and baryon mass fractions. For a given halo mass, escape 
fractions can vary up to a factor of two, depending on the initial setup of the idealized halo. In 
a cosmological setting, we find that the time-averaged photon escape fraction always exceeds 25% 
and reaches up to 80% in halos with masses above 10 8 solar masses with a top-heavy IMF. The 
instantaneous escape fraction can vary up to an order of magnitude in a few million years and tend 
to be positively correlated with star formation rate. We find that the mean of the star formation 
efficiency times ionizing photon escape fraction, averaged over all atomic cooling (T vlr > 8000 K) 
galaxies, ranges from 0.02 for a normal IMF to 0.03 for a top-heavy IMF, whereas smaller, molecular 
cooling galaxies in minihalos do not make a significant contribution to reionizing the universe due 
to a much lower star formation efficiency. These results provide the physical basis for cosmological 
reionization by stellar sources, predominately atomic cooling dwarf galaxies. 

Subject headings: cosmology: theory — galaxies: formation — galaxies: dwarf — galaxies: high- 
redshift — radiative transfer 



1. INTRODUCTION 

In most calculations of cosmological reionization, it has 
to be assumed that the product of star formation effi- 
ciency and hydrogen ionizing photon e scape fraction is 
> 0.01 ( e.g. lGnedinll200a |Cenll2003allbL IWvithe fc Loebl 



2003allbT IVenkatesan et all 120031: ISomerville et all 120031: 
Chiu et all 120031: iHaiman fc Holder! 120031: ICiardi et^al I 
2001 iSokasian et all ffifll l200i IWvithe fc Cenl [20071 : 

Srbinovskv fc Wvithel l2008') in order to reionize the uni- 
verse early enough to be consistent with the Wilkin- 
son Microwave Anis o tropy Probe (WMAP ) observations 
(jSpergel et all l2007t iKomatsu etall l2008f). if sta r s pro - 
duce the majority of ionizing photons. lGnedin| (2000) 
suggested that, based on Local Group dwarf galaxies, the 
star formation efficiency at high re dshift is ~ 4%, con- 
sistent with theoretical work s (e.g. iKrumholz fc McKed 
120051: IKrumholz fc Tanll2007fh Thus, unless star forma- 
tion efficiency is unusually high (i.e., > 10%) at high 
redshift, this requirement seems to suggest a high ion- 
izing photon escape fraction from high redshift galaxies, 
/esc > 10%, may be necessary in the context of stellar 
reionization within the standard cold dark matter model. 

However, / csc > 10% is by no means the norm, at 
least from observations at lower redshifts. Approxi- 
matel y 6% of ionizing radiation escape from the Milky 
Way (Bland-Hawt horn fc Malonevl Il999t iPutman et al.l 
120031 ). For local starburst galaxies (Hurwitz et all (|1997l ) 
gave / esc < 0.032,0.052,0.11 (2a) for Mrk 496, Mrk 
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1267 . and IRAS 08339 +6517 (< 0.57 in the case of Mrk 
66). iDeharveng et alJ (2001) gave an escape frac tion of 
/esc < 0.06 2 for Mrk 54. iHeckman et alJ (|2001h found 
/esc < 0.06. iBergvall et al.l (|2006f) found /U. ~ 0.04-0.1 
for a local extrem e starburst dwarf, the Blue Compact 
Galaxy Haro 11. IChen et all (|2007h placed a 95% up- 
per limit of /esc < 0.075 for star forming regions hosting 
gamma-ray b ursts at z > 2 . For Lyman break galax- 
ies at z ~ 3 IShaplev et all (|2006h found / csc ~ 0.03. 
ISiana et all (|2007l ) examined starburst galaxies at z ~ 
1.3 and found that less than 20% have relative escape 
fractions 4 near unity and a global absolute escape frac- 
tion of <0.04. However in this study, so me galaxies 
have u pper limits of / eS c,rei as low as 0.08. Ilnoue et alJ 
(2006) recently compiled available observations over a 
wide range of redshift and come to the conclusion that 
there is a trend of increasing / csc with redshift from 
< 0.01 at z < 1 to ~ 0.1 at z > 4. 

Some theoretical works also seem to point to low val- 
ues for / esc . Theoretical models of our own Galaxy with 
realistic star formation mode by iDove et "ahl (I2000T) give 
an estimate of / esc ~ 0.06 or lower. ICiardi et al.l (2002) 
studied the effect of gas inhomogeneities on / esc for Milky 
Way like galaxies using simulations with radiative trans- 
fer and concluded that / OS c depends on the density struc- 
ture as well as star formation rate (SFR). In addition, 
porosity in the interstellar medium (ISM) that is caused 
by SN mechanical feedback may provide additional chan- 
nels in which radiat ion could escape, thu s increasing the 
UV escape fraction (jClarke fc Oe vll2002T ) . iWood fc Loebl 

4 This quantity is the absolute escape fraction that is corrected 
for dust attenuation at 1500 A and is a common measure in obser- 
vations of escaping Lyman continuum. 
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(2000 ) argued that f e sc < 0.01 for galaxies at z ~ 10. 
iRicotti fc Shul (|2000h found that / osc > 0.1 only for 
halos of total mass less than ~ 10 7 M© at z > 6 with 
/ esc d ropping precipitously for larger halos. iFuiita et all 
(2003), using ZEUS-3D simulations, found that / esc < 
0.1 from disks of d warf starburst galaxies of tota l mass 
M > 1O 8-1O M . iRazoumov fc Sommer-Larsenl (|2006l 
120071 hereafter RS06 and RS07) concluded that / esc 
= 0.01-0.1 in several young galaxies with M v j r = 
10 12 ~ 13 M© at z ~ 3 in their simulations with radiative 
transfer using adaptive ray tracing, agreei n g with the ob- 
servat ions presented in llnoue et al.l (|2006t ) . iGnedin et al.l 
(|2rM hereafter GKC08), using detailed hydrodynamic 
simulations with radiative transfer, found / csc ~ 0.01 — 
0.03 for galaxies of total mass M > 10 11 M at z = 3 — 5. 

In the relevant redshift range for cosmological reion- 
ization, z = 6 — 15, most of the ionizing radiation due 
to stars comes from dwarf galaxies of M < 10 8 ~ 9 M w in 
the s tandard cold dark matter model (iBarkana fc Loebl 
2001). This purpose of this paper is to investigate how 
/csc depends on various physical parameters expected 
in realistic cosmological settings using detailed adaptive 
mesh refinement (AMR) simulations coupled with adap- 
tive 3D ray-tracing radiative transfer, focusing on galax- 
ies with total mass M = 10 6 5 - 1O 95 M at z > 6. We 
study the dependence of / esc on four physical parame- 
ters: mass of the galaxy, spin parameter, baryonic mass 
fraction and tu rbulent energy. Ou r work is an extension 
of GKC08 and lFujita et al.l (j2003h . by exploring the de- 
pendence of /csc on some important physical parameters 
aforementioned, by including lower mass galaxies and by 
employing very high resolution (O.lpc) simulations. In 
addition, importantly, we allow for disparate lifetimes of 
stars of different masses, of their subsequent explosive 
energy (for those stars that become supernovae) input 
into the ISM. 

We first describe our radiation hydrodynamics simu- 
lations of isolated and cosmological halos in There 
we also describe our algorithm for star formation and 
feedback that considers a multi-phase ISM and resolved 
molecular clouds. In £ 13.11 and §3.2[ we present the star 
formation rates and history and the resulting escape frac- 
tion of UV radiation, respectively. Next in $4l we further 
compare our results with previous work and discuss ef- 
fects from any neglected physical processes in our simu- 
lations and the implications of our results on reionization 
scenarios. We summarize our work in the last section. 

2. RADIATION HYDRO DYNAMICAL 
SIMULATIONS 

We use the Eulerian ada ptive mesh refinement (AMR) 
hydrodynamic code Enzo (jBrvan fc Normanlll997i 119991 : 
lO'Shea et ail l2004h to investigate the escape of ioniz- 
ing radiation from early dwarf galaxies. We have mod- 
ified the code to include a time-dependent implemen- 
tation of radiative t ransfer that adaptively ray traces 
from point sources (lAbel fc WandehH 120021 lAbel et all 
120071 1 Wise fc Abel |20Q8a| ). Enzo uses an n-bod y adap- 
tive particle-mesh solver (|Efstathiou et al.lll985l ) to fol- 
low the dark matter (DM) dynamics. It solves the 
hydrodynamical equation using t he second-order accu- 
rate piecewise parabolic method ([Woodward fc Colellal 
119841: iBrvan et alj Il995f ). while a Riemann solver en- 
sures accurate shock capturing with minimal viscosity. 



We use the nine-species (H I, H II, He I, He II, He III, 
e — , H2 , H2 + , H~ ) non- e quilibrium chemistry model in 
Enzo (|Abel et al.l Il997t lAnninos et all 119971) and the 
H 2 cooling rates from iGalli fc Pallal (|1998t ). We adopt 
the cosmological parameters from th e third year "mean" 
WMAP results (jSpergel et al.ll2007l ) of fl M = 1 - ft a = 
0.24, fl b h 2 = 0.0229, and h = 0.73, where Q A , fl M , and 
fib are the fractions of mass-energy contained in vacuum 
energy, cold dark matter, and baryons, respectively. Here 
h is the Hubble constant in units of 100 km s _1 Mpc -1 . 
For the cosmological simulations, we use a § = 0.76 and n 
= 0.96, where as is the rms of density fluctuations inside 
a sphere of radius 8ft. -1 Mpc and n is the scalar spectral 
index. 

We first describe our simulation setup of our suite of 
idealized dwarf galaxies. Next we detail our cosmologi- 
cal simulations that we use to check the validity of our 
assumptions in the idealized cases. Last we explain our 
star formation and feedback model. 

2.1. Isolated Halo Setup 

In order to quantify the dependence of / csc on dwarf 
galaxy properties, we run a total of 75 idealized dwarf 
galaxy simulations, in which we vary the total halo mass 
M v ; r , spin parameter 



(1) 



baryon mass fraction fa = Mg as /M v i r , and gas turbu- 
lence /turb = Vrms/V c . Here L, E, and M gas are the 
total angular momentum, energy, and gas mass in the 
halo. When initializing the halos, E is the total grav- 
itational energy of the halo, and during the analysis, 
E = M«j ms /2. v rms is the three-dimensional rms gas 
velocity, and 

v. -(^) V2 

\ ^ vir / 

is the circular velocity of the halo, where r V j r is the virial 
radius. We define the virial radius and thus mass as the 
radius of a sphere that encloses an overdensity of 200 
times the cosmic mean density. 

We calibrate our choice of such parameters by the fol- 
lowing previous work. 

1. Spin parameter — It has b e en well-established 
dBarnes fc Efstathioul Il987t lEisenstein fc Loebl 
1995) that the distribution of the halo spin pa- 
rameter is log-norm al with a mean of .042 and 
deviation of 0.5 (e.g. iBullock et al.ll2001at ). 

2. Baryon fraction — Radiative feedback from the 
first stars and subsequent episodes of star for- 
mation generate outflows that reduce the baryon 
fraction up to a factor of 3 below the cos- 
mic aver age (Flb/flM) in halos with M v - lv 



10 8 Mp (IRicotti fc Gnedin) 1 20051: iWise fc Abe 
2008aT IRicotti et alj 120081: iMesigner fc Diikstn 
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3. Turbulent energy — Turbulence increases during 
the assembly of a cosmological halo through both 
virialization and mergers. Typical Mach numbers 
in halos with an adiabatic state of equation is ~0.25 
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TABLE 1 

COSMOLOGICAL HALO PROPERTIES 



High Luminosity Normal Luminosity 



# 


M vir 


fb 


f~tuvb 


A 


-^grid.O 


A^part 


{/"esc } 


<SFR) 




max(SFR) 


{fese } 


(SFR) 




max(SFR) 




[M e ] 












[Mq yr" 


l ] 


[Mq yr" 1 ] 




[Mq yr" 




[Mq yr" 1 ] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 




(10) 


(11) 


(12) 




(13) 


1 


1.6 x 10 7 


0.11 


0.56 


0.057 


118 3 


5,144 


0.38 


3.7 x 10" 


-5 


1.7 x 10~ 4 


0.065 


1.5 x 10- 


-4 


4.1 x 10~ 3 


2 


2.6 x 10 7 


0.12 


0.48 


0.038 


130 3 


8,482 


0.77 


3.0 x 10- 


-4 


1.1 x 10" 3 


0.13 


6.7 x 10- 


-4 


4.8 x 10" 3 


3 


2.7 x 10 7 


0.13 


0.61 


0.061 


132 3 


8,809 


0.42 


1.0 x 10' 


-3 


3.7 x 10" 3 


0.11 


9.5 x 10" 


-4 


3.1 x 10" 3 


4 


5.0 x 10 7 


0.13 


0.57 


0.013 


154 3 


16,103 


0.67 


1.1 x 10" 


-3 


4.1 x 10~ 3 


0.44 


4.1 x 10- 


-3 


1.6 x 10~ 2 


5 


1.7 x 10 s 


0.14 


0.67 


0.044 


240 3 


52,831 


0.78 


2.6 x 10" 


-2 


7.6 x 10" 2 


0.38 


1.6 x 10- 


-2 


6.5 x 10" 2 


6 


3.8 x 10 s 


0.12 


0.60 


0.027 


80 3 


1,222 


0.44 


2.3 x 10" 


-2 


8.9 x 10~ 2 


0.39 


5.0 x 10- 


-2 


1.1 x lO" 1 


7 


6.7 x 10 s 


0.12 


0.58 


0.052 


84 3 


2,184 


0.72 


5.3 x 10" 


-2 


2.2 x lO" 1 


0.36 


5.4 x 10" 


-2 


1.7 x 10- 1 


8 


8.7 x 10 8 


0.13 


0.54 


0.025 


90 3 


2,783 


0.47 


4.8 x 10" 


-2 


1.9 x lO" 1 


0.31 


6.3 x 10- 


-2 


2.2 x lO" 1 


9 


1.2 x 10 9 


0.13 


0.62 


0.049 


98 3 


3,863 


0.51 


5.7 x 10" 


-2 


1.9 x lO" 1 


0.40 


5.3 x 10- 


-2 


1.1 x lO" 1 


10 


4.0 x 10 9 


0.13 


0.70 


0.069 


156 3 


12,796 


0.62 


1.7 x 10" 


-1 


5.1 x lO" 1 


0.54 


2.2 x 10- 


-1 


4.2 x lO" 1 



Note. — Low and high luminosity models have 2,600 and 26,000 ionizing photons per stellar baryon. Col. (1): Halo number. 
Col. (2): Virial mass. Col. (3): Baryon mass fraction. Col. (4): Turbulent energy. Col. (5): Spin parameter. Col. (6): 
Number of grid cells on the top-level AMR grid. Col. (7): Number of dark matter particles inside the virial radius. Col. (8,11): 
Time-averaged UV escape fraction. Col. (9,12): Time-averaged star formation rate. Col. (10,13): Maximum star formation rate. 



(iNorman fe Brvanl 119991 : iNagai fe Kravtsovl 120031 
iWise fe Abell l2007aF ~ In halos with gas that can 
efficiently cool, the halo cannot no longer virialize 
through gas heating, becoming more turbulent to 
virialize. The combination of lower temperatures 
and greater turbulent motions leads to turbulent 
Mach numbers above unity and clos e to the circu- 
lar velocity of the halo, i .e. /t U rb = 1 (jWise fe Abell 
[2007al iGreif et al.ll2008h . 

We hence fix two of the following parameters (A, fb, /turb) 
= (0.04, 0.1, 0.25) and change the remaining parameter 
and halo mass in each calculation. We allow these halo 
parameters to vary as: 

• logio Af vir = 6.5, 7, 7.5, 8, 8.5, 9, 9.5; 

• A = 0, 0.02, 0.04, 0.06, 0.1; 

• h = 0.05, 0.075, 0.1, 0.15; 

• /tmb = 0, 0.25, 0.5, 1. 

For example, when varying /;,, we model halos with all 
of the quoted virial masses and fb while fixing A and 
/turb to 0.04 and 0.25, respectively, resulting in a total 
of 28 calculations. We repeat this method for varying 
the spin parameter and turbulent energy. We did not 
evaluate a M v i r = 10 6 5 M Q halo in two cases, fb = 0.05 
and fturb = 1, because they did not condense to form 
stars. 

The simulation box has a side length of 10 times the 
virial radius at a virialization redshift z v i r = 8 and a 
top grid resolution of 128 3 . We refine the grid struc- 
ture if the gas density becomes 1.5 times greater than 
the mean gas density times a factor of 2 l , where I is 
the AMR refinement level. We also refine to resolve 
the local Jeans length by at least 4 cells. Cells are 
refined to a maximum AMR level (ranging from 9-12) 
that results in a spatial resolution of ~0.1 proper par- 
sec, which is required to model the for mation of the 
D-type front at small scal es correctly (|Whalen et al.l 
l2004t iKitavama et all [20041) . Here we sample the ini- 
tial Stromgren sphere by several cells across, which has 



a radius of ~1 pc if we assume an ionizing luminosity of 
10 50 photons s" 1 and a density of 1800 cm" 3 (see 35T3]) . 

We perform the calculations in comoving coordinates. 
The gravitational potential is computed from an NFW 
density profile, 

Pnfw{t) = j—, — \7~T~\ — / — v>> ( 3 ) 
(r/r s )(l + r/r s y 

where r s is a characteristic inner radius, and p s is the 
corresponding inner density this is chosen so the total 
DM mass within r v ; r is (1 — /f>)A/ V ir- The core radius can 
be expressed in terms of the concentration parameter c = 
r s /r v ii- We set c = 6, wh ich is in concordance with the fit 
of iBullock et al . (2001b). In addition, we consider self- 
gravity from the time-dependent baryon density field. 

The grid boundaries are periodic, which has no effect 
on the evolution of the galaxy because our simulation 
box is much larger than the virial radius. The initial 
configuration of each halo is an isothermal sphere, i.e. 
p cx r~ 2 , with a constant density core with radius r s . 
The initial temperature is set to the virial temperature. 
Rotational energy is modelled as a solid body rotator, 
and turbulent energy is initialized as random velocities 
sampled from a Maxwellian distribution with a temper- 
ature T = (/turb^c) 2 /im p /fc, where [i is the mean molec- 
ular weight in units of a proton mass m p . In most of our 
models, turbulent energy is dom i nant oy er rotational en- 
ergy, as shown in [Wise fc Abell (|2007a[ ). We start our 
calculations at redshift 8. The surface of the halo is in 
pressure equilibrium with the IGM that has a density of 
QbPa where p c = 3Hq/8itG is the critical density of the 
universe. 

2.2. Cosmological Halo Setup 

In reality, these halos reside in a cosmological environ- 
ment, where adjacent filamentary structures may inhibit 
the escape of ionizing radiation. To investigate any ef- 
fects of a cosmological setting, we run two sets of AMR 
cosmological simulations with periodic boundary condi- 
tions and 384 3 particles and grid cells. The first set of 
simulations (named M6.5-8) is set up to study halos with 
masses between 10 6 5 and 10 8 M Q with a side length of 
1.75 comoving Mpc, resulting in a DM particle mass of 
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2720 Mq. The second set of simulations (named M8.5- 
9.5) has a comoving side length of 8 Mpc and focuses on 
halos with masses between 10 85 and 10 95 M Q , resulting 
in a DM particle mass of 2.72 x 10 5 M Q . Our choice of 
resolution and box size allows halos in the quoted mass 
range to be resolved by at least 1000 DM particles. We 
initialize simulations M6.5-8 and M8 .5-9.5 at a redshift 
68 and 51, respectively, with graf ic ( Bci tschi ngerll200"l[) 
and different random phases. AMR refinement criteria 
are based on the same quantities as the idealized halos, 
now also include DM density, but with an overdensity 
criterion of 4 instead of 1.5 for both gas and DM. We al- 
low the grid structure to be refined in the entire volume. 
These simulations arc stopped at redshift 8, which is the 
same redshift the idealized halos are initialized. Note 
that GKC08 found no redshift dependence in / esc . At 
the final redshift, simulations M6.5-8 and M8.5-9.5 have 
132,456 and 114,602 AMR grids and 5.63x 10 8 (826 3 ) and 
3.11 x 10 s (677 3 ) unique computational cells, respectively. 

We use an adiabatic equation of state in these cosmo- 
logical simulations because the inclusion of atomic line 
and molecular hydrogen cooling would create star form- 
ing, cold, dense molecular clouds. Currently it is not 
feasible to treat >, 300 point sources with our ray trac- 
ing scheme; therefore we suppress any gas cooling and 
thus star formation in these large-scale simulations and 
take the following approach. 

At the final redshift for 10 selected halos with at least 
1000 DM particles, we extract a sub-volume with its 
AMR hierarchy intact and a box length of 10 times its 
virial radius. The halo masses range from 1.6 x 10 7 to 
4.0 x 10 9 Af Q , and they do not undergo a major merger 
in the 100 Myr after z = 8. We take the resolution of 
a level 3 AMR grid in the large-scale simulation to be 
the top-level resolution of the sub-volume, ranging from 
80 3 to 240 3 . Additional details of these halos are listed 
in Table [TJ Now each sub- volume is its own simulation, 
and we consider the nine-species non-equilibrium chem- 
istry model and star formation and feedback. We only 
refine the grid structure within the inner quarter of the 
box that is centered on the halo of interest. We calcu- 
late the gravitational potential with isolated boundary 
conditions, and the potential is assumed to be zero at 
the boundaries. We use inflowing and periodic boundary 
conditions for baryons and DM particles, respectively. 
We restrict star formation to within 1.5 virial radii of 
the halo center. We adjust the baryon and DM veloci- 
ties so that the mass-averaged velocity of the simulation 
is zero. We also set the fractional abundances of (H II, 
He II, He III, H~, H 2 , H+) = (1.2 x 10~ 5 , 10" 14 , 10~ 17 , 
2 x 10~ 9 , 2 x 10~ 6 , 3 x 10~ 14 ). Although these val- 
ues are not the equilibrium values, the non-equilibrium 
chemistry solver quickly converges to the correct values 
during the first few timesteps of the calculation. 

These halos may represent a lower limit of / csc because 
previous episodes of star formation may have decreased 
the baryon mass fraction and the column densities of ad- 
jacent filaments; thus, / esc may be larger in the realistic 
case that examines both Population II and III star for- 
mation during the entire halo assembly history. However, 
we can use the results from our idealized halos to under- 
stand how / osc behaves with different galactic properties 
and make the according adjustments to our estimates of 



/ csc in these extracted cosmological halos. 

2.3. Star Formation and Feedback 

We model star formation throu gh a modified version o f 
the extensively used algorithm of lCen fc Ostrikerl (|1992h , 
where a stellar particle is created when all of the follow- 
ing criteria are met: (i) an overdensity of 10 7 , equivalent 
to 1800/x -1 cm -3 at redshift 8; (ii) a converging gas ve- 
locity field (V • v < 0); (iii) rapidly cooling gas, i.e. the 
cooling time t coo \ is less than the dynamical time ta yn = 
l/(Gp) 1 / 2 . Since we ensure that the local Jeans length is 
always sampled by at least 4 cells, a single cell will never 
be Jeans unstable 5 , thus we do not consider the Jeans 
unstable criterion outlined in the original work. When a 
cell satisfies the above criteria, we define the star forming 
cloud as a sphere with td yn = 3 Myr (p c i — 1000^ cm" 3 ) 
and radius i? c /, where a fraction c* = 0.07/ C oid of the gas 
is gradually converted into a stellar particle, whose final 
mass = c Jr (47r/3)p c /i? 3 ; . Here / co id is the fraction of 
gas with T < 10 4 K within the sphere. This treatment of 
cold gas accretion is simila r to the multi-phase model o f 
star formation specified in lSpringel fc Hernquistl {2003). 
We take this approach because we resolve the molecular 
cloud instead of employing a subgrid model. The stellar 
mass in the particle increases as 



dM _ 
~dT 
or equivalently 



yn 



t-ti 



td 



yn 



exp 



-(t-*i) 



td 



yn 



M{t) (t-U) + t dyn 
— — = 1 - — exp 

ITl* ^dyn 



-(t-tj) 
tdyn 



(4) 



(5) 



We then replace the sphere with a uniform density p c i — 
(l-c*)/(G4J. 

Our cho ice of the star formation ef ficiency is in agree- 
ment with lKrumholz fc McKei (|2005f) . who find that this 
fraction in a free-fall time is 0.077a~°- 68 A^" - 32 , where 
a v i r is a virialization parameter and is near unity, and 
M is the turbulent Mach number. We justify our choice 
of a threshold idyn with inferred dynamical times from 
observations of star forming regions are approximately 
700 kyr, a nd star formation occurs for several dynamical 
times fe.g. lTan et aLll2006ft . 

In order to minimize the number of radiative sources 
in the calculations, we enforce a minimum star particle 
mass of 10 _5 M gas , where M gas is the gas mass inside the 
halo, and merge stellar particles if they are within i? me rgo 
pc. To keep the number of sources below ~100, we have 



found the relation R n 



10(M vir /10 8 M Q ) a5 pc works 



well. The merged particle inherits the center of mass as 
its position, mass-averaged velocity, and total mass of its 
two progenitors. We found that the UV escape fraction 
is not significantly dependent on this parameter and only 
varies by A(f csc )/f csc — 0.1 when we change R mcTgc by 
a factor of 0.1, 0.5, and 2. 

Radiative stellar f eedback is modelled wit h an adaptive 
ray tracing scheme (|Abel fc Wandeltl l2002t ) that is cou- 
pled to the hydrodynamics, energy, and chemistry solver 

5 Refinement based on Jeans length usually occurs in moderate 
refinement levels (i.e. 5-8). With primordial gas cooling, dense 
gas does not cool below 200 K, and at a density of 1800 cm -3 , 
the corresponding Jeans length is ~5 pc, which is well above our 
resolution limit. 
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of Enzo (jWise fc A bel 20081, Wise et al., m prepara- 
tion) . It is parallelized with MPI and runs on distributed 
and shared memory machines. Each stellar particle is 
a radiative point source and has a constant luminos- 
ity of A 7 = 26,000 ionizing photons per stellar baryon 
over its lifetime. This ionizing luminosity is appropri- 
ate for a Salpeter initial mass function (IMF) with a 
metallicity of 4 x lO~' l Z p ) and lower and upper mass cut- 
offs of 1 and 150 M Q (|Schaeredl200l . Similarly high 
specific luminosities may also be applicable in gas with 
Z >, 10~ 25 Zq at high-redshift whose cooling is limited 
to the CMB temperature. This raises the Jeans mass 
of such molecular clouds and may result in a top-heavy 
IMF with a typica l stellar mass of tens of solar masses 
( Smith et al.l [2008). To compare our work to lFuiita et al.l 
|2003|) and GKC08 in the "control halos" of the ideal- 
ized setup and cosmological halos, we use iV 7 = 2,600 
that is suitable for a solar- metallicity Salpeter- like IMF. 
The photons are evenly distributed among 192 rays, i.e. 
HEALPix level 2, emanating from its associated radia- 
tion source. As the rays propagate from the source, the 
solid angle that is sampled by each ray increases as r 2 . 
When this solid angle becomes larger than 20% of the 
cell area that it is tra nsversing, the ray splits accord- 
ing to the formalism of I Abel Sz Wandelti The rays cast 
in these simulations have an energy of 21.5 eV, which 
is th e mean ionizin g photon energy from our fiducial 
IMF (Schacrcr 2003j). We have experimented with pho- 
ton means energy between 14.0 and 30.0 eV and found 
that our results are not sensitive to this energy. Aver- 
aged over 10 Myr, this results in an ionizing luminosity of 
10 46 and 10 47 erg s _1 M Q _1 for our chosen stellar lumi- 
nosities. We consider only hydrogen photo-ionization in 
these simulations. We do not consider dust absorption 
as GKC08 showed it is only important redward of the 
Lyman break. We model the H2 dissociating radiation 
between 11.2 and 13.6 eV with an optically-thin 1/r 2 ra- 
diati on field c entered on each source with luminosities 
from iSchaererl . We use the H2 photo-dissociation rate 
coefficient fc^iss = 1.1 x 10 8 Flw s _1 , where Flw is th e 
H 2 flux in units of erg s" 1 cm~ 2 Hz -1 (jAbel et al.lll997t ). 

After the stellar particle has lived for 4 Myr, we model 
supernova (SN) feedback by injecting (dM/dt)(At) 6.8 x 
10 48 erg Mq 1 of thermal energy at every timestep of the 
finest AMR level into a sphere with a diameter of 1 pc. 
If the grid cell that the star particle resides in is larger 
than 1 pc, we inject the energy into the surrounding 27 
grid cells for numerical stability. For SN feedback, we 
evaluate dM/dt at time, t - 4 Myr, to model the SN from 
the stars that were "born" 4 Myr prior to the current 
time. We eliminate the star particle when it has an age 
of 10 Myr, corresponding to the lifetime of OB-stars that 
produce the vast majority of ionizing photons. 

2.4. UV Escape Fraction 

We measure the fraction of escaping UV radiation by 
comparing the number of photons that travel beyond the 
virial radius to the current stellar luminosity. In other 
words, the integral form of the UV escape fraction is 

/ esc (i) = -1— f F vv (t)dS, (6) 

where L\jy and F\jy are the UV stellar luminosity and 
UV flux, respectively, and <S is the spherical surface at 




log Halo Mass [ivy 

Fig. 1. — The mean of the time-averaged star formation rates 
(SFR) for idealized halos. The thin and thick lines correspond to 
runs with a normal (7V 7 = 2, 600) and top-heavy (7V 7 = 26, 000) 
IMF, respectively. The error bars show the minimum and maxi- 
mum values at each halo mass. The average SFR for cosmological 
halos with a top-heavy (normal) IMF are represented with dia- 
monds (plus signs). The dotted lines show the SFR for a constant 
star formation efficiency, /* = M*/M gas , for /j, = 0.1 and lasting 
for 100 Myr. 

T = r v i r . The equivalent quantity that is calculated in 
the simulations is 

, , v = N i( r vir <r < rvir + c dtp) 

Usc[) - Y,Ni{r<cdt p ) ' [<) 

where iV» is the number of photons contained in the i-th 
photon package, r is the total distance that the pack- 
age has travelled, and c dt p is the distance travelled by 
a photon within a timestep dt p . The numerator sums 
the number of photons that have travelled farther than a 
virial radius within the last timestep, and the denomina- 
tor is the number of photons emitted in the last timestep, 
which can be easily counted when the photon packages 
are generated from the point source. It is possible to 
compare these two quantities that are temporally sepa- 
rated by the light-travel time (t <~ 3 kyr) , which is much 
smaller than the dynamical times of the star formation 
regions, to which the change in SFR is directly related 
(see Eq. HH5]). Also notice that a photon at a distance 
r v i r from its source does not necessarily correspond to 
the surface S; however this is a good approximation to 
/esc because the stars are centrally concentrated in its 
halo galaxies. 

3. RESULTS 

3.1. Star Formation and Feedback 

Here we present our results of our simulations that con- 
sider the self-consistent treatment of star formation and 
feedback in high-redshift dwarf galaxies with masses be- 
tween 10 6 5 and 10 9 ' 5 solar masses. We first discuss the 
amount and trends with halo mass of the star formation 
in these galaxies. Then we illustrate the effects of radia- 
tive feedback on the gas morphology and self-regulation 
of star formation. Last we describe the star formation 
histories of such galaxies. 

3.1.1. Star Formation Rates 

The escape fraction of UV radiation from halo is pri- 
marily dependent on the strength of the intrinsic ionizing 
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Fig. 2. — The maximum SFR for idealized halos with a normal 
(thin dashed) and top-heavy IMF (thick dashed) and cosmological 
halos with a normal (crosses) and top-heavy IMF (diamonds). 





luminosity and how it propagates through the gas into 
the IGM. The halos undergo an initial collapse and ensu- 
ing starburst. Afterwards star formation is self-regulated 
by radiative and supernovae feedback. We find that most 
of the UV radiation escapes during the most intense 
episodes of star formation, thus an important quantity 
to investigate is the maximum SFR of the halo, as well 
as the average. We plot the time-averaged and maximum 
SFR of the idealized and cosmological halos with a top- 
heavy and normal IMF in Figures Q] and [21 respectively. 
We also plot SFRs if a constant star formation efficiency, 
/* = M*/Mg as , is assumed for a baryon mass fraction of 
0.1 and a starburst lasting for 100 Myr. 

The maximum SFRs are anywhere between 2-10 times 
the average SFR. The error bars show variances up 
to an order of magnitude in the star formation for a 
given halo mass in the idealized runs. For halos with 
M vir > 10 7 M Q , SFRs steadily increase as SFR oc M v - lr a 
with a ~ 2. Idealized halos with M vlr = 10 6 5 and 10 7 
Mq have similar SFRs because they experience only one 
episode of star formation that blows away most of the 
gas in the halo, quenching any further star formation. In 
a more realistic cosmological setup, these outflows only 
delay later star formation because the gas reservoir can 
be reple nished through gas ac cretion from filaments and 
mergers (|Wise fc Abelll2008a| ). This accretion causes the 
SFRs to be higher up to a factor of 100 in cosmological 
halos up to Myir = 10 8 ' 5 Mq, which can be seen in Fig- 
ures [T] and [2j Above this halo mass, SFRs interestingly 
level off at ~0.1 Mq yr _1 , corresponding to an average 
star formation efficiency /* = 0.05 — 0.1. 

With a normal IMF, average and maximum SFRs in 
cosmological halos with M v ; r < 10 8 M Q increase by an or- 
der of magnitude. In these halos, outflows generated by 
radiative feedback is the main force in regulating SFRs. 
Therefore one expects less self-regulation of star forma- 
tion, and thus higher SFRs, with a normal IMF when 
compared to a top-heavy one. In more massive halos, 
only outflows generated by SN explosions can escape 
from the potential well. Hence decreasing the specific 
stellar luminosity does not significantly affect the gas 
mass inside the halo, leading to similar SFRs in the top- 
heavy and normal IMF runs. 




Fig. 3. — Gas density projections, weighted by density squared 
for idealized halos. The field of view is 2r v ; r . The colormap spans 
5 orders of magnitude from 10 -3 to 10 2 cm -3 . 



3.1.2. Stellar Feedback 

In the dwarf galaxies that we consider, stellar feed- 
back from H II regions and SN explosions creates galac- 
tic outflows that expel significant amounts of gas from 
the halo. Our simulations capture the formation of a 
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Fig. 4. — Gas density projections, weighted by density squared for cosmological halos. Halos 1—10 have logio(A^vir) 
8.2, 8.6, 8.8, 8.9, 9.1, and 9.6. The field of view is 2r v i r . The color scaling is the same as Figure |31 



7.2, 7.4, 7.4, 7.7, 



D-type ionizatio n front that accelerates matter up to 
35 km s" 1 fe.g iFranco et al.lll99d IWhalen et aD l2004t 
iKitavama etal1l2004h lAbel et al.ll2007fr which is greater 
than the escape velocities of halos with M v ; r < 1O 85 M0 
at redshift 8. Above this mass, SNe feedback, which gen- 
erates outflows up to ^300 km s _1 , becomes the dom- 
inant force in expelling gas. Feedback strongly affects 
the gas structures within the dwarf galaxies, which we 
illustrate in density-squared projections of gas density of 
idealized and cosmological halos in Figures [3] and 2] at 
10, 50, and 100 Myr after the initial starburst. 

The transition from complete gas "blow-out" to "blow- 
away" is apparent in the idealized halos (Fig. [3]) at 
M> 1O 8J5 M . In halos with M vil = 10 6 - 5 ~ 7 Mq, the D- 
type front is unbound from the halo and expels most of 
the gas with it. Only a small fraction returns to the halo 
center after 100 Myr. In halos with M vir = 10 7 - 5_8 M Q , 
the blow-out becomes less effective at completely evacu- 
ating the halo because of the deeper potential well. As 



stated before when the D-type front velocity is less than 
the escape velocity of the halo, radiative feedback cannot 
create unbound galactic outflows. Nevertheless, it still 
plays a key role in shaping the small-scale (/ <, 100 pc) 
variations in the ISM density and temperatures. 

The cosmological halos provide a more realistic en- 
vironment and halo setup — inflows from filaments, mi- 
nor mergers, turbulent velocity fields, an evolving DM 
potential — in studying the effects of star formation and 
feedback. Despite the more realistic gas inflows in such 
halos, we find a similar transition from blow-out to blow- 
away at M vir = 10 8 - 5 M©. The halos with M vir < 10 8 M Q 
are not evacuated like their idealized counterparts be- 
cause of additional gas accretion from adjacent filamen- 
tary structures and minor halo mergers 6 . In response, 

6 In even less massive halos (M v ; r ~ 10 6 Mq), Pop III radiative 
feedback can completely evacuate its host halo in a cosmological 
environment lAbel et a l. 2007; Yoshida et al. 2007; Johnso n et al.l 
[20071; IWIse fc AbeW2008ah . 
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star formation is self-regulating, in which the stellar radi- 
ation reduces the amount of cool and dense gas available 
for star formation. In the least massive halos, star for- 
mation ceases for tens of millions of years as the relic 
H II region cools and re-condensed. We further discuss 
star formation rates in i J3.1.3ll3"T.4l 

Above 10 s Mq, radiative feedback does not destroy 
all of the density enhancements, in particular one cen- 
tred in the halo center, and stars continue forming in 
the presence of nearby star clusters. The infalling sub- 
halo (M vir ~ 10 7 M Q ) in Halo 6 at t = 100 Myr provides 
an interesting contrast in the halo mass dependence of 
radiative feedback, where the main halo is centrally con- 
centrated while the subhalo is being blown apart by ra- 
diative feedback. In these more massive halos, star for- 
mation is still self-regulated but not suppressed as there 
is a persistent cold gas reservoir available for star forma- 
tion, as seen in the evolution of gas density in the right 
column of Figure |U 

3.1.3. Star Formation History of Idealized Halos 

We plot the star formation histories of the idealized 
control halos (f b = 0.1, A = 0.04, / turb = 0.25) in the 
top panels of Figure [SJ Star formation in halo masses 
below 10 8 M Q is episodic. This occurs when radiation 
driven stellar outflows expel most of the gas from the 
halo, quenching star formation. Once the gas is no longer 
irradiated, it can cool and fallback into the halo center, 
forming stars once again. This occurs up to 5 times in 
these halos in the 100 Myr, more often in more massive 
halos. 

When halo masses increase to between 10 7 5 and 10 85 
Mq, star formation continues in the presence of galactic 
outflows; however it soon halts once the cool gas reser- 
voir in the halo is depleted. In Figure [5j one sees in halos 
with these masses experience at least one period of qui- 



escence, where a substantial fraction of gas cannot form 
star-forming molecular clouds because it is contained in 
outflows. Similar to the lower-mass halos, star formation 
recommences after gas falls back into the halo. 

Above 10 9 M Q , the bursting behavior of star forma- 
tion is entirely eliminated. The potential well of the 
halo has become deep enough (V c > 35 km s _1 ) to 
contain any outflows created by over-pressurized H II 
regions. The transition circular velocity from baryon 
"blow-out" to "blow-away" is similar to the velocities 
generated in D-type ionization fronts. Furthermore, un- 
like the lower mass halos, the roughly constant SFRs 
are approximately the same between the high- and low- 
luminosity models at 5 x 10~ 2 and 10 _1 Mq yr _1 for 
M vir = 10 9 and 1O 9 - 5 M , respectively. 

Thin disk formation is absent in halos that experience 
baryon blowout because the outward motions created in 
H II regions disrupt any global organized rotation. In 
halos with M v - U > 1O 9 M0, the majority of gas is re- 
tained within the halo, which then settles into a thin 
disk. We caution that is probably not strong evidence 
for a transition from a globally turbulent ISM to disk 
formation because of our idealized setup of solid-body 
rotation. We can better investigate any morphological 
trends in cosmological halos that are described in the 
following section. 

3.1.4. Star Formation History of Cosmological Halos 

We show the star formation histories for the 10 se- 
lected cosmological halos in the top panels of Figure El 
The SFR in cosmological halos is less episodic than the 
idealized halos because the cosmological gas accretion 
replenishes the galaxy, balancing the effect of any blow- 
away. Nevertheless star formation is still self-regulated, 
which can be seen in the SFR being suppressed up to 
an order of magnitude after a local maximum in SFR. 
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Same as Figure[5]but for cosmological halos. The halo numbers are labeled in parentheses. 



Afterwards the SFR recovers again because of the re- 
duced radiative feedback. This modulation occurs in all 
cosmological halos studied here. 

In lower mass halos with M v ; r < 1O 8 M with a top- 
heavy IMF, star formation is still completely eliminated 
for 10-20 Myr. This does not occur in the same ha- 
los with a normal IMF because the galactic outflows are 
weaker, making it easier for the cosmological accretion to 
penetrate through the outflows. Hence star formation is 
continuous in these cases. Similar to the M v ; r > lO 9 Af 
idealized halos, SFRs are not significantly affected by 
a lower specific stellar luminosity in cosmological halos 
with M v j r > 10 8 A/q. There are some deviations by a 
factor of a few in the evolution of the SFR in these ha- 
los, but the time-averaged rates are largely unaffected. 
One clear example of this behavior occurs in Halo 6 
(M vir = 3.8 x 10 8 M Q ) at t = 70 Myr. In the top-heavy 
IMF case, radiative feedback from the initial burst sub- 
dues star formation from t = 50-80 Myr by a factor of 
5. In contrast, the feedback from a normal IMF does 
not suppress star formation, which continues at a rate of 
~1O _1 M0 yr _1 . Furthermore, galactic disk formation, 
which is seen in isolated halos with M v - lr > 10 9 M Q , is 
absent in all of the cosmological cases. 



3.1.5. Comparison to Previous Work 

The SFRs in the M vir = 1O 95 M halos are 
^10 _1 Mq yr _1 and are similar to the least massive 
(M vir - 10 10 M Q ) halos in GKC08. However for a direct 
comparison of SFRs and escape fractions, we must cither 
use the simulations with iV 7 — 2,600 or consider the SFRs 
from to p-heavy IMF simul ations as effectively 10 times 
greater. iFuiita et al.l (|2003f ) characterizes star formation 
with a star formation efficiency = M+fMd, where Md 
is the total gas disk mass. If we equ ate the total gas 
mass in our halos to the disk mass of IFuiita et all our 
M V ir = 10 9 M Q idealized halos are comparable to their 
high-redshift M d = 1O 8 M models with /* = 0.06. 

3.2. UV Escape Fraction 

A fraction of the stellar radiation emitted from the 
stars described in the previous section escape into the 
IGM. This quantity is an essential para meter in bo th 
semi-analytical reionization models (e.g. iCenl l2003aH bl; 
iWvithe fc Loebll2003at lHaiman fc Holderl [20031) a nd N- 
body reionization simulations (e.g. Ciardi e talj [20031 ; 
ISokasian et aT1[2^IIh^et~aT] [2006: Tr ac k Cen|[200l . 
First we discuss the trends of / csc with halo mass of the 
idealized halos. We then describe the behavior of / esc 
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Fig. 7. — Escape fraction of ionizing photons from idealized 
halos with a top-heavy IMF (solid line) and normal IMF (circles) 
and cosmological halos with a top-heavy IMF (filled diamonds) and 
normal IMF (open diamonds) as a function of halo mass. 



with respect to important physical parameters: turbu- 
lent energy, halo spin parameter, and baryon mass frac- 
tion. Next we discuss the results of / esc in halos that 
have been extracted from cosmological simulations. We 
end this section discussing anisotropic H II regions and 
its relevance in the UV escape fraction. 

3.2.1. Idealized halos 

We plot the evolution of / csc of the idealized halos in 
the bottom panels of Figure [5] Before presenting our re- 
sults of /esc in these halos, we would like to warn that one 
should not take the values of / esc as absolute because of 
the lack of gas accretion and adjacent filamentary struc- 
tures. Nevertheless it is appropriate to study the tempo- 
ral variations and changes with physical parameters. 

In halos with M v i r < lO 7 Af0, / csc behaves in a pre- 
dictable manner, in which star formation is confined 
within a small region and happens in bursts. Thus most 
of the escaping radiation experiences similar gaseous 
structures when propagating outwards. In the bursts, 
the / csc quickly increases when the H II region breaks 
out from the halo, which can take up 10 Myr after the 
starburst starts. The escape fraction peaks when the 
SFR is at its maximum and varies anywhere from 10 -3 
to unity, depending on the column density of the sur- 
rounding medium. 

In more massive halos, the escape fraction is highly de- 
pendent on preceding star formation and its effect on the 
gas structure. For example in halos with M vlr = 10' 5 
- 10 8 5 Mq, SFRs that differ by an order of magnitude 
can produce similar values of / CS c. In the initial star- 
burst, the radiation needs to escape through a relatively 
smooth r~ 2 density profile and / 0S c > 0.5. Then condi- 
tions quickly erode from this smooth initial state to a tur- 
bulent multi-phase ISM because of both radiative cooling 
and stellar feedback. Now radiation must escape through 
a clumpy and turbulent ISM, whose porosity can increas e 
/csc for a given luminosity or SFR (|Clarke fc Oevi r2002). 

Generally only radiation from the strongest starbursts 
can escape, whereas during suppressed periods of star 
formation very little (/ es c < 0.01) radiation escapes. 
As SFRs become more consistent in halos with M v ; r > 
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Fig. 8. — Escape fraction of ionizing photons as a function of 
maximum star formation rate from all simulated halos, showing 
a strong dependence for a given halo mass. Idealized halos with 
equal halo masses arc plotted with the same symbols. Idealized 
halos with a normal IMF are plotted as large circles ranging from 
10 6 5 to 1O 9 ' 5 M0 from left to right. Cosmological halos with a top- 
heavy and normal IMF are plotted as filled and open diamonds, 
respectively. 



1O 9 M , the escape fraction varies between 0.1 and 1.0 
with less variation than the lower mass halos. In addi- 
tion, we find that radiation from stars that form on the 
outskirts of the thin disk have a higher escape fraction 
than stars that are deeply embedded within the galaxy, 
which is similar to the results of GKC08. 

In Figure we plot the average UV escape fraction 
as a function of halo mass, where the error bars indicate 
the standard deviation of / csc over the 11 models for 
each halo mass. Halos with M v ; r < 10 s M Q have / csc ~ 
0.5 ± 0.2. For a normal IMF in these low-mass halos, 
/ esc drops to values below 0.1. Escape fractions steadily 
increase with halo mass to ~0.8 above this mass for both 
top-heavy and normal IMFs. 

As seen in Figure / OS c is highest when the SFR is 
at its maximum, max(SFR). To depict this strong de- 
pendence on the strength of the starburst, we show / esc 
for every simulated halo as a function of max(SFR) in 
Figure [8j grouped by halo mass. For instance in the 
M vir = 10 7 - 5 M Q halos, / CS c increases from 0.2 to 0.8 
when max(SFR) is higher by an order of magnitude. We 
have also plotted the idealized halos with a normal IMF 
as large open circles. 

3.2.2. Dependence on halo parameters 

Although the absolute values of / esc in idealized ha- 
los are inaccurate for reasons discussed in the previ- 
ous section, we can utilize the trends with physical halo 
parameters — baryon mass fraction, spin parameter, and 
turbulent energy — to estimate / esc in cosmological halos 
other than the ones s tudied in this paper. For example, 
IWise fc Abell (|2008ah found that the first stars lowered 
the baryon mass fraction to < 0.1, but the halos from 
our adiabatic cosmological simulations have /& ~ 0.13. 
We can then use the trends in / esc with respect to fb to 
estimate the effects of neglected processes on the escape 
fraction. We show these trends for all idealized halos in 
Figure \9\ 
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Fig. 9. — UV escape fraction versus halo mass plotted for varying values of (a) turbulent energy, (b) spin parameter, and (c) baryon 
mass fraction. 



Turbulent energy — Initial turbulent energy should be 
unimportant in the more massive halos as radiative and 
SN feedback co ntinuously alters the ISM a nd further stirs 
turbulence (e.g. IScalo fc Elmegre cn 2004) , perhaps over- 
powering any pre-existing turbulence. There are no con- 
sistent trends with turbulent energy between 10 - 10 9 5 
M G . However in the M vir = 1O 7 5 M halos, the f tur b = 
case has / esc = 0.8 that is up to twice the escape frac- 
tion of the other halos. This occurs because turbulent 
flows are not present to disrupt the initial halo collapse, 
resulting in a stronger initial starburst and escalating 
the average escape fraction. The fturb — case with 
M v ir = 10 8 M© is only marginally higher than the other 
halos by ^25%. This anomaly is not present in other 
halo masses. 

Spin parameter — Similar to the lack of apparent 
trends with turbulent energy, spin parameter has lit- 
tle effect on the UV escape fraction. One exception is 
the non-spinning case above 10 7 M©, having the largest 
escape fraction. This arises from a strong central col- 
lapse without being rotationally supporting, resulting 
in a stronger initial starburst. However this scenario 
is improbable in r eality because of abundan t cosmolog- 
ical tidal torques <|Hoyldll949t lPeebleslll969h . culminat- 
ing in a log-normal spin para meter distribution that has 
a mean value of 0.0 42 (e.g. H arncs fc Efstathioul H987t 
iBullock et afll2lioTan . 

Baryon mass fraction — Out of the three physical pa- 
rameters studied, the baryon mass fraction plays the 
largest role in regulating star formation and the escaping 
radiation. In halos with M V - 1T < 1O 8 M0, the fb = 0.15 
case has / eS c lower than the other gas fractions by a factor 
of 2. The increased neutral column density that the ra- 
diation must transverse before breaking out of the halo 
causes the decrease. Interestingly above 10 8 M Q , this 
trend reverses, and higher gas fractions result in larger 
escape fractions. Here the gas-rich halos can sustain a 
higher SFR, and the corresponding higher ionizing lumi- 
nosity can escape the halo more freely than in the gas- 
poor (fb < 0.1) halos, even though the radiation has a 
larger column to ionize. 



3.2.3. Cosmological halos 

The star-forming cosmological halos represent a more 
realistic setup of high-redshift, starburst dwarf galaxies, 
necessary to calculate a robust value of the UV escape 
fraction. In the lower panels of Figure we show the 
escape fraction of the cosmological halos as a function 
of time. During the first starburst that originates from 
the initial halo collapse, instantaneous / esc values widely 
vary from less than 10~ 3 up to nearly unity. Most im- 
portantly the escape fraction depends on the current star 
formation and any of its time variations. As seen in Fig- 
ure [5] the variations in SFR corresponds to variations 
in /csc on similar timescales, on order of a dynamical 
timescale of a molecular cloud, ~l-3 Myr. In halos 2 
and 4, a constant SFR exists for 30 Myr in the top-heavy 
IMF case. The cumulative luminosity from this stellar 
cluster gradually ionizes the IGM and photo-evaporates 
most clumpy material. Eventually the escape fraction 
approaches unity after ~10 Myr of irradiation. 

As with the idealized halos, the escape fraction is de- 
pendent on the previous star formation history, which 
can photo-evaporate dense, clumpy material in the halo. 
Imagine two equal-mass stellar clusters, forming in the 
same neighborhood. The one that forms first will pre- 
ionize the ISM in some fraction of solid angle, allowing 
the radiation from the second cluster to escape into the 
IGM more easily and further raising the escape fraction. 

We compare / csc from cosmological halos with their 
idealized counterparts in Figures [7] and [8] with respect to 
virial mass and maximum SFR, respectively. The halos 
that host a top-heavy IMF have escape fractions within 
the scatter of the idealized halo, even though adjacent fil- 
amentary structures can absorb most of the radiation in 
their line of sight. Their / csc values range from 0.4 to 0.8. 
We also find that / esc at half and twice the virial radius 
are within 10% of the value at r V j r , similar to GKC08, 
because the nearby clumps and adjacent filaments do not 
contribute a significant solid angle for absorption. 

A normal IMF lowers the escape fraction by A(/ csc ) = 
0.05-0.4. It is interesting that f csc is approximately 0.4 
with a normal IMF over 2 orders of magnitude in halo 
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20.51 21.23 21.96 22.69 20.37 21.10 21.83 22.56 -2.00 -1.50 -1.00 -0.50 21.28 21.98 22.67 23.36 15.96 18.38 20.80 23.23 -2.00 -1.50 -1.00 -0.50 



f ^9^9 vw 




20.52 21.05 21.59 22.12 15.71 17.79 19.86 21.93 -2.00 -1.50 -1.00 -0.50 21.82 22.32 22.82 23.32 17.02 19.07 21.12 23.17 -2.00 -1.50 -1.00 -0.50 
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20.22 20.67 21.12 21.57 20.00 20.48 20.97 21.45 -2.00 -1.50 -1.00 -0.50 20.98 21.82 22.68 23.50 15.79 18.31 20.84 23.37 -2.00 -1.50 -1.00 -0.50 



20.56 20.84 21.11 21.39 16.01 17.84 19.28 20.91 -2.00 -1.50 -1.00 -0.50 20.75 21.68 22.60 23.53 20.41 21.41 22.41 23.41 -2.00 -1.50 -1.00 -0.50 
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Fig. 10. — Full-sky maps of total column density (first and fourth columns), H I column density (second and fifth columns), and mass- 
weighted neutral hydrogen fraction (third and sixth columns) 50 Myr after the initial starburst in the cosmological halos. Halos 1—10 have 
log 10 (Af v i r ) = 7.2, 7.4, 7.4, 7.7, 8.2, 8.6, 8.8, 8.9, 9.1, and 9.6. The observer is at the center of mass of the halo, and the quantities shown 
out to the virial radius. The colormaps for column density are unique for each halo; however the electron fraction colormaps have a fixed 
range of 0.01—0.3. Halos 2 and 4 are almost completely ionized, and halos 6, 7, and 10 represent clear examples of anisotropic H II regions, 
which escape from the halo where the column density is lowest. 



mass and maximum SFR, starting with M vir = 1O 7 - 5 M . 
Halos below this mass threshold (i.e. efficient atomic 
cooling) are affected considerably by a normal IMF, 
where / esc drops by a factor of a few to values of 0.05-0.1. 

3.2.4. Anisotropic H II Regions 

In our radiation hydrodynamics simulations, radiative 
feedback greatly affects the gas dynamics inside the ha- 
los. Any UV radiation will first escape from the halo in 
the direction with the least H I column density. This 
creates anisotropic H II regions with the radiation pref- 
erentially escaping through these channels. Any adja- 
cent filamentary str uctures, which provide the halo with 
a cold, dense flow (iNagai et al.l [20031: 
iDekel fc Birnboiml 12009: iWise fc Abel 



Keres et al.l 120051: 
20073), and cold 



and clumpy ISM are the major components in absorbing 
outgoing radiation. 

We create full-sky maps of total and H I column density 
and mass-averaged neutral fractions of the cosmological 
halos with a top-heavy IMF 50 Myr after the initial star- 
burst in Figure ITD1 These maps are created by casting 
adaptive rays (|Abel &: W andelt 2002) from the center of 
mass of the halo to the virial radius. Each ray tracks 
the total and neutral column density along its path. We 



then reconstruct a full-sky map of these quantities us- 
ing the spatial information contained in the HEALPix 
formalism. The ionized regions (blue) match well with 
the areas with the lowest column density (black and dark 
blue). The transitions from neutral to ionized in these 
maps are abrupt, depicting how any neutral blobs can 
completely absorb the radiation in its solid angle. How- 
ever, in the solid angles where the ionized fraction is 
close to unity, almost all of the radiation can escape, i.e. 
/ csc = 1 in those directions. The simulations of RS06, 
RS07, and GKC08 also see the majority of solid angles 
having / esc equal to either unity or zero. 

At t — 50 Myr, halos 1 and 3 are in a quiescent state, 
and the previously ionized ISM and IGM is recombining. 
The relic H II region are still visible in the full-sky maps. 
Here the transition from ionized and neutral regions are 
gradual because of the longer recombination times in dif- 
fuse gas. Halos 2 and 4 have been almost completely 
ionized at the end of a starburst. Notice the cometary 
tails created by photo-ionization of a clumpy m edium 
(e.g. lSusa fc Umemura|[200llWhalen et al.ll2008ft at the 
northern section in halo 2 and the middle-bottom region 
in halo 4. In the former case, the neutral shadow, point- 
ing southwest, of the small overdensities are clear in the 
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neutral fraction map. In halos 5-10, stellar radiation es- 
capes through a large fraction (> 25%) of solid angle, 
which correspond to the path of least neutral column 
density. 

3.2.5. Comparison to Previous Work 

Clearly these escape fractions contradict the results 
from GKC08, where their lowest mass halos with M v ; r ~ 
10 10 M Q have / csc = 1(T 5 - 10~ 2 at z = 3 - 5. From the 
examples given there, the stars are born within a galac- 
tic disk, unlike our simulations, which could be causing 
their smaller escape fractions. In 21 we further discuss 
possible causes of this discrepancy. Although our galax- 
ies are smaller by a mass factor of >100 than the ones 
simulated in RS06 and RS07, it is worth reiterating that 
they find / csc = 0.01 — 0.1 at redshift 3. Both groups 
also find that / e sc can significantly vary from galaxy to 
galaxy and in diffe rent lines o f sight to each galaxy. Es- 
cape fractions in iFuiita et a l. (2003) exceed 0.2 in their 
high-redshift dwarf galaxies if the starbursts are intense, 
i.e. their /* > 0.06 models, similar to our / esc values in 
the most massive halos. 



4. DISCUSSION 

The fraction of ionizing radiation that escapes high- 
redshift galaxies is an important value to quantify in 
order to characterize cosmological reionization. In this 
section, we first discuss the differences and similarities, 
along with their respective causes, between our results 
and previous work. Next we discuss any implications 
and agreements with semi-analytic reionization models. 
We last elaborate on possible influences from neglected 
physical processes. 

4.1. Possible Physical and Computational 
Dependencies 

Our results suggest that / OS c is higher than 0.1 in such 
galaxies, but we must first understand why our / esc val- 
ues differ from the more massive galaxies presented in 
the s imula tions of iRazqumov fe Sommer-Larsenl (|2006l . 
12001 and iGnedin et all (120081) . Below we discuss five 
possible causes. 

Galaxy morphology and environment — The galaxies 
simulated here exhibit an irregular morphology and did 
not form a rotationally-supported disk, whereas the sim- 
ulations of RS06, RS07, a nd GKC08 a l l stud ied disk 
galaxies. As mentioned in IFuiita et al.l (|2003l ). a tur- 
bulent and clumpy ISM may allow for radiation to es- 
cape more easily into the IGM. Compared to a disk 
configuration, there are more low-density escape routes, 
i.e. porosity, for t he radiation becaus e of the clumpy 
nature of the gas (|Clarke fc Oevl [20021 ). Furthermore, 
the stellar orbits are not confined within a disk and can 
reach the outskirts of the galaxy at apocenter, where gas 
densities are much lower than the galactic center and 
ionizing radiation can escape into the IGM at a much 
greater rate. Perhaps this irregular morphology is only 
initially present in dwarf galaxies because radiative feed- 
back sh ifts the angular mom entum distribution to higher 
values ([Wise fc Ab el 2008a). A fraction of this material 
will return to the galaxy and aid in disk formation in 
higher mass galaxies, such as the ones studied in RS06, 



RS07, and GKC08. Furthermore in rare halos, the inter- 
secting filaments tend to be dense and thin, where a halo 
with equal mas s at lower redshift w ould be contained in a 
large filament (|Ocvirk et al.ll2008t ). This redshift depen- 
dency on halo environment could affect / C sc at r > r V j r , 
but should not be apparent at r v ; r . 

Galactic mass — GKC08 find that / osc increases as 
much as 3 orders of magnitude from a halo mass of 10 10 to 
lO 11 -/!//©. The cold H I disk in their lower mass galax- 
ies tend to be more vertically extended, which brings 
about this precipitous drop. If this trend continues to 
lower masses, the values of / csc would be negligible in 
our galaxies, but we see no such continuation. Our re- 
sults suggest that the radiative feedback in dwarf galaxies 
with V c < 35 km s" 1 have a profound impact on the ISM 
by driving outflows and preventing any disk formation, 
increasing / esc . Perhaps this circular velocity is a critical 
turning point below which / osc is large (> 0.1) because of 
dynamical effects from radiative feedback. In more mas- 
sive galaxies, SN feedback provides the main impetus for 
driving outflows. 

Star formation rates — One can argue that our simula- 
tions only capture the initial collapse and an artificially 
strong starburst in high-redshift galaxies; however, after 
~20 Myr of star formation, this collapse is reversed by 
stellar feedback and SFRs stabilize. In addition our most 
massive galaxy (M v ; r = 4 x 10 9 M Q ) has an average SFR 
of 0.2 Mq yr _1 , which is comparable to the SFRs found 
in the M vir - 10 10 M Q galaxies in GKC08. The star 
formation efficiencies (M*/M gas ) of our galaxies range 
from 5-10%, also agreeing with previous galaxy forma- 
tion simulations. Hence we do not think our treatment of 
star formation affects the SFRs and resulting UV escape 
fraction. 

Stellar IMF — We experimented with both normal 
{N 1 = 2, 600) and top-heavy {N 1 = 26, 000) IMFs and 
found that a normal IMF causes / esc values to be lower 
by 10-75%. The differences in how the gas and ensuing 
star formation reacts to radiative feedback most likely 
causes this spread. The normal IMF is similar to the 
ones used in all previously quoted theoretical works, and 
thus we do not expect our choice in IMF to cause our 
high UV escape fractions. One minor shortcoming of 
our simulations is that we keep the SN feedback strength 
equal in our normal and top-heavy IMFs. Although this 
helps us localize the cause of any changes to a different 
specific luminosity, the strength of SN feedback should 
also affect the gas distribution, SFR, and / csc . 

Resolution — We can afford a spatial resolution of ~0.1 
physical pc in all of our simulations. However in larger 
galaxies, it becomes increasingly difficult to maintain this 
resolution, especially in cosmological simulations. The 
highest resolution simulation in RS06 and RS07 use a 
gravitational softening length of 330 h~ x pc, and their ra- 
diative transfer grid is refined so that 10 particles occupy 
each cell. GKC08 use AMR si mulations wi t h a m aximum 
resolution of 50 physical pc. iFujita et al.l (|2003l ) have a 
fixed resolution of 0.28 pc in their z = 8,M d = 1O 8 M 
simulation, similar t o our resolut ion. The / esc values of 
our simulations and IFuiita et al.l are in agreement. Reg- 
ular and giant molecular clouds have ra dii of 2-20 and 
10-60 pc (see iMac Low fc Klessenl I2004L for a review), 
respectively, and are usually accounted with subgrid 
physics in cosmological simulations. However our sim- 
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ulations capture the turbulent and clumpy nature of the 
ISM, while resolving such molecular clouds. As discussed 
before, the turbulent ISM and its po rosity may play a ke y 
role in allowing radiation to escape (|Clarke fc Oevll2002h . 
In addition, high-resolution simulations with radiative 
tran sfer can model the fragmentation of ionizatio n fronts 
(e.g. lFuiita eTal]|2003t IWhalen fc Normanll2008[ ). which 
can further assist in boosting UV escape fractions. 

4.2. Impact on reionization models 

In reionization models, whether it be based on Prcss- 
Schetcher formalism or N-body simulations, the prod- 
uct /osc = /* x /esc ultimately dictates the absolute 
number of photons that escape from each halo into the 
IGM. When calibrating these models against WMAP 
observations (jSpergel et all l2007t iKomatsu et all l2008h 
and Lyq transmi ssion f rom z ~ 6 quasars (e. ; 
iBolton fc Haehnelt] [20071 : iSrbinovskv fc Wvithd l200l 
the values of > 0.01 are usually required. Recently, 
ISrbinovskv fc Wvithd found that / C sc must lie within the 
range 0.1-0.25 with a best-fit of /* = 0.11 if all halos 
with masses > 10 9 M Q contribute to reionization. Fur- 
thermore, they conclude that even the smallest atomic 
line cooling halos (M vir > 10 8 M Q ) have / osc ~ 0.05. Our 
results support this scenario where low-luminosity galax- 
ies are the main contributors to reionization. 

We plot /* S c fr° m idealized and cosmological halos in 
Figure [TTJ We calculate /* using the initial gas mass 
in halo at the initial time; thus we are overestimating 
/* up to a factor of ~1.5 because the halos experience 
mass accretion in the 100 Myr that we have simulated. 
Recall that none of the halos presented here undergo a 
major merger. In our least massive cosmological halos 
with M v i r < 3 x 10 7 M Q , /* sc increases rapidly with re- 
spect to halo mass from 10~ 3 to 10~ 2 because atomic 
hydrogen line cooling becomes efficient in these halos. 
/* sc is also an order of magnitude higher than the ide- 
alized cases. As discussed in .1.11 the additional gas 
accretion from filaments results in higher SFRs or equiv- 
alently /*. In atomic cooling halos (T v ; r > 8000 K), we 
find /* sc to always be > 0.01, sufficiently high enough 
to meet the "critical" value to match reionization con- 
straints, with an average of 0.033 and 0.026 for top-heavy 
and normal IMFs, respectively. When integrated over a 
luminosity function with a faint-end slope of -1.7 (e.g. 
iBouwens et al.|[2007t ). the average of /* sc in atomic cool- 
ing halos is 0.029 and 0.021 for top-heavy and normal 
IMFs, respectively. If we include all halos, this average 
drops to 4.1 x 10~ 3 because the low-mass halos are more 
abundant and cannot cool and form stars as efficiently. 

In Figure [T2l we show the behavior of /* sc with halo 
mass in idealized halos grouped by halo parameters, sim- 
ilar to Figure [9] There is no apparent trends with respect 
to turbulent energy; however, halos with higher spin pa- 
rameters result in smaller values of /* sc in high mass halos 
that form a rotationally supported disk. The halos with 
fb = 0.1 have the highest values of /* sc , and gas-poor 
halos are consistently lower in all halo masses. 

4.3. Effects from the Ultraviolet Background 

Our simulations accurately track the evolution of the 
H II regions and how they breakout into the IGM. Here 
we discuss an important process that was neglected in our 
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Fig. 11. — The product of stellar mass fraction (M s tar /-Mgas ) 
and escape fraction of ionizing photons (/esc) from idealized ha- 
los and cosmological halos with a top-heavy IMF (filled diamonds) 
and normal IMF (open diamonds) as a function of halo mass. Ide- 
alized halos with a normal IMF are plotted as large circles. This 
product is a key quantity in semi-analytical reionization models 
in determining the amount of escaping radiation per collapsed gas 
fraction. 



simulations, external feedback from the UV background 
(UVB). 

Our models sample halos that are primarily dependent 
on H2 cooling (T V - 1T < 8000 K) and atomic line cooling 
in more massive halos. Gas condensation in the lower 
mass H2 cooling halos can be delayed by a H2 disso- 
ciating UV ( Lyman- Werner) background between 11.2 
and 13.6 eV dMachacek e t alJl2001H Wise fc Abe]| l2007bt 
lO'Shea fc Normanl |200^ T The Lyman- Werner back- 
ground thus increases cooling times in the centers of such 
halos. As a result, the minimum mass of a star-forming 
halo increases with the Lyman- Werner background in- 
tensity. This will not affect the global halo properties 
but may suppress the SFR in these halos. The Lyman- 
Werner background becomes less of an issue in atomic 
line cooling halos as Lya cooling provides ample amounts 
of free electrons for H2 co oling, and they bec ome self- 
shielding to this radiation ()Wise fc Abelll2008aD . 

We now consider any photo-heating from a hydro- 
gen ionizing UVB, which can partially suppress gas ac- 
cretion and thu s star formation in dwarf galaxies (e.g. 
EfstathioulH99l lShaDirolt~aLlfl99l fThoul fc Weinberg 



1996). Here the Jeans "filtering mass" ([Gnedin fc Hui 
19981 ) accurately describes the minimum halo mass that 



can cool and collapse given an IG M thermal history 
(|Gnedinl l2000t iWise fc Abel l2008bl) . The Jeans filter- 
ing mass is calculated by analyzing the linear evolution 
of overdensities in the presence of Jeans smoothing. Be- 
cause the low-density IGM has a dynamical time on the 
order of a Hubble time, it slowly reacts to any photo- 
heating, and the filtering mass can be thought of some 
running time-average of the Jeans mass. Thus at some 
rcdshift, t he Jeans and filtering ma ss may differ sub- 
stantially. iThoul fc Weinberg! (jl996) originally showed 
photo-heating could totally suppress any star formation 
in low-redshift dwarf galaxies with U c <,35 km s~ x and 
somewhat lower SFRs in galaxies up to 100 km s _1 . 
At z >3, halos are less susceptible to negative feedback 
from photo-heating because of their higher average densi- 
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Fig. 12. — Relative changes in the product of stellar mass fraction (M atar /Af gas ) and escape fraction of ionizing photons (/ eS c) from all 
idealized models grouped by initial halo parameters: turbulent energy (left), spin parameter (middle), and baryon mass fraction (right) 
when compared to the "control halos" with (/ tur b; A, /;,) = (0.25, 0.04, 0.1). 



ties and thus cooling rates (jDijkstra et alJl20"04[ ). Closer 
inspection of collapsing halos in radiation hydrodynam- 
ics simulations also shows high-redshift, low-mass halos 
can collapse and is in deed regulated by the filtering mass. 
Prior to reionization, Gnedin| (|2000D showed that the fil- 
tering mass smoothly increases from ~10 6 before any 
star formation occurs to 10 9 Af© at z = 6. IWise k Abel 
(|2008bf ) calculated the filtering mass in regions that were 
ionized by Population III stars and found that it gradu- 
ally increases the filtering mass to ~3 x 10 7 M Q at z ~ 15 

around biase d regions. 

Recently iMesigner fe Diikstral (|2008f) studied how 
photo-heating from an inhomogeneous UVB suppresses 
gas cooling and star formation. They found that at 
z = 10 an UVB intensity of J21 ~ 0.1 is necessary to 
completely suppress gas collapse in M V i r = 10 8 M Q ha- 
los, where J21 is in units of 10 -21 erg s _1 cm~ 2 Hz -1 
sr _1 . In their semi-numerical simulations at z = 10, 
they showed this intensity occurs in a volume fraction 
of < 1%, suggesting that star formation is widespread 
in these low-mass halos even in the advanced stages of 
reionization. They also illustrate how the collapsed gas 
fraction steadily decreases with increasing UVB. Even 
at z — 7 in a 10 8 M Q halo, approximately 40% of the 
gas is retained when compared to the no feedback case. 
This gas mass loss will directly affect star formation 
and perhaps t he subsequent UV es cape fraction. Using 
the results of [Mcsigncr & Diikstra] and the dependence 
of /* sc on gas fraction (see Fig. [12j) , our results from 
cosmological halos can be adjusted to account for ex- 
ternal feedback from the UVB. For example, at z = 10 
with their fiducial model, their semi-numerical simula- 
tions show that J21 ~ 0.01 is most common, which de- 
creases the gas fraction of 10 8 , 10 85 , 10 9 Af Q halos to 
0.2, 0.6, 0.85, respectively, of the gas fraction without 
any UVB. From Fig. [T2l we can see that /* sc drops by 
^50% when the gas fraction is <0.075 in idealized halos. 
We would expect a similar decrease in halos that lose gas 
from UVB feedback, resulting in /* sc ~ 0.02 for halos in 
this mass range. 



4.4. Prior Star Formation Episodes 

We assumed that the halos were unaffected by any 
prior stellar feedback, but in this paper, we have also 
shown that low-mass halos are susceptible to gas "blow- 
out" caused by radiative feedback. This effect should 
be e vident starting with Population III star formation 
(e.g. lYoshida et atl 120071 : IWise fc Abel l2008a| ) in halos 
with ci rcular velociti es up to V c ~ 35 km s _1 . For ex- 
ample, IWise fc Abell found that gas fractions of dwarf 
galaxies with M v i r ~ 3 x 10 7 M Q at redshift 15 have 
been decreased to < 0.1. The cosmological halos pre- 
sented here, which have /& ~ 0.13, have / csc ^0.5 with 
a top-heavy IMF; however if fb were lower from previ- 
ous gas "blow-out" , the trends shown in Figure [9] sug- 
gest that the escape fraction in these halos should be 
even higher, but their values of /* sc (Fig. [T2f will de- 
crease because of lower SFRs. To further understand and 
quantify the consequences of these neglected processes, 
semi-analytical reionization models, semi-numerical or 
N-body + radiative transfer simulations are needed, as 
discussed in §4.21 In future work, we plan to study this 
in detail by following the hierarchical assembly of high- 
redshift dwarf galaxies in cosmological simulations, in- 
cluding both metal-free and metal-enriched star forma- 
tion and feedback. 

5. SUMMARY 

We presented results from an extensive suite of very 
high resolution (0.1 pc) AMR radiation hydrodynamics 
simulations that focus on the UV escape fraction from 
isolated halos and cosmological halos. The latter cases 
were extracted from two large-scale cosmological simula- 
tions. These simulations accurately track the evolution of 
H II regions while including radiative and SNe feedback 
on the surrounding medium. The halo masses studied 
in this work span from 3 x 10 6 to 3 x 10 9 M Q . We used 
the isolated halos to gauge if the escape fraction depends 
on the turbulent fractional energy, spin parameter, and 
baryon mass fraction of the halo. The cosmological ha- 
los provide a good estimate of realistic escape fractions 
of high-redshift dwarf galaxies. We also investigated how 



16 



WISE & CEN 



a top-heavy IMF and a normal IMF affects the escape 
fraction. Our main findings in this paper are 

1. Radiative feedback from massive stars, primarily 
arising from D-type fronts, is most effective at ejecting 
gas from halos with masses less than 10 8 - 5 M Q (V c = 
35 km s" 1 ). 

2. Radiation preferentially escapes through channels 
with low column densities, creating anisotropic H II and a 
highly varying / osc along different lines of sight, agreeing 
with previous work. 

3. For a given halo mass, escape fractions in isolated 
halos have a standard deviation of 0.2, arising from dif- 
fering physical halo parameters. The largest effect comes 
from the baryon mass fraction, where / osc is lower in gas- 
rich halos with masses below 10 8 M . In gas-rich halos 
with masses larger than 10 8 M Q , star formation is more 
intense, overcoming any large neutral fraction it must 
ionize and resulting in higher values of / osc than halos 
with smaller gas fractions. 

4. High-redshift dwarf galaxies with M v j r > 1O 7 M , 
a top-heavy IMF, and irregular morphology have 0.25 < 
/osc < 0.8, which we determined from our simulations 
of cosmological halos. A normal IMF decreases / csc to 
0.05-0.1 in halos with M vir < 1O 7 5 M and / esc ~ 0.4 in 
more massive halos, which have maximum SFRs span- 
ning almost 2 orders of magnitude. 

5. Escape fractions are dependent not only on the cur- 
rent SFR but on the photo-heating and dispersion of gas, 
following feedback from previous episodes of star forma- 
tion. Values of / osc can vary up to an order of magnitude 
in a few million years, i.e. the dynamical time of a molec- 
ular cloud on which variations in SFRs can occur. 



6. The mean of the product of star formation efficiency 
and ionizing photon escape fraction, averaged over all 
atomic cooling (T v ; r > 8000 K) galaxies, ranges from 0.02 
for a normal IMF to 0.03 for a top-heavy IMF. Smaller, 
molecular cooling galaxies in minihalos are not significant 
contributors to reionizing the universe primarily because 
of a much lower star formation efficiency in minihalos 
than in atomic cooling halos. 

The high escape fraction of UV radiation has im- 
portant implications on reionization, allowing a large 
amount of radiation from low-luminosity dwarf galaxies 
to freely propagate into the IGM. Although our simula- 
tions miss the entire assembly of halos and their complete 
star formation history, they are robust in following the 
dynamics of star formation and feedback during the 100 
Myr studied and suggest that escape fraction in these 
galaxies are larger than previously assumed. 
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